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ABSTRACT 


We present chemical implications arising from spectral models fit to the Her- 
schelfiHFl spectral survey toward the Orion Kleinmann-Low nebula (Orion KL). We 
focus our discussion on the 8 complex organics detected within the HIFI survey uti¬ 
lizing a novel technique to identify those molecules emitting in the hottest gas. In par¬ 
ticular, we find the complex nitrogen bearing species CH 3 CN, C 2 H 3 CN, C 2 H 5 CN, and 
NH 2 CHO systematically trace hotter gas than the oxygen bearing organics CH 3 OH, 
C 9 H 5 OH, CH 3 OCH 3 , and CH 3 OCHO, which do not contain nitrogen. If these com¬ 
plex species form predominantly on grain surfaces, this may indicate N-bearing organ¬ 
ics are more difficult to remove from grain surfaces than O-bearing species. Another 
possibility is that hot (7’ km ~ 300 K) gas phase chemistry naturally produces higher 
complex cyanide abundances while suppressing the formation of O-bearing complex 
organics. We compare our derived rotation temperatures and molecular abundances to 
chemical models, which include gas-phase and grain surface pathways. Abundances 
for a majority of the detected complex organics can be reproduced over timescales 
> 10 5 years, with several species being under predicted by less than 3cr. Derived 
rotation temperatures for most organics, furthermore, agree reasonably well with the 
predicted temperatures at peak abundance. We also find that sulfur bearing molecules 
which also contain oxygen (i.e. SO, S0 2 , and OCS) tend to probe the hottest gas to¬ 
ward Orion KL indicating the formation pathways for these species are most efficient 
at high temperatures. 

Subject headings: astrochemistry — ISM: abundances — ISM: individual objects (Orion KL) 


— ISM: molecules 
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1. Introduction 


Hot cores, the birth sites of massive stars, are the most prodigious emittere of molecular line 


radiation in the Galaxy, harboring high abundances of organic molecules (Herbst & van Dishoeck 


20091) . Many of the organics detected in these regions are composed of 6 atoms or more, which are 


often labelled “complex” by astronomers. Though the chemistry in these regions remains poorly 
understood, theoretical studies argue that the evaporation of icy grain mantles, precipitated by the 


formation of protostars, pi 

ays a key role in the 

roduction of many complex organic mo 

ecules 

(see e.g. 

Millar et al. 

1991 

Chamlev et al. 

1992; 

Caselli et al. 

1993; 

Charnlev et al. 

1995 

1997 

Rodgers & Chamlev 

2001 ; 

Garrod & Herbst 

2006; 

Garrod et al. 

200 

8 ). However, the extent to 


which these species form on grain surfaces as opposed to in the gas phase remains unclear. 

The Orion Kleinmann-Low nebula (Orion KL) is a high mass star forming region 
that provides a superb template for constraining the chemical processes active within these 
environments. Orion KL harbors one of the most organic rich hot cores in the Milky Way 


(B lak e et al. 


1987). Moreover, the so-called compact ridge, a 


clumps residing ~ 12" south-west of the hot core (Wang et al. 


group of externa 


2011 


Favre et al. 


y heated dense 


201 lal) . is equally 


rich in organic species. Relative to one another, the hot core and compactjidge are more abundant 


in organics that contain nitrogen and oxygen, respectively (IBlake et al. 


2005 


Guelin et al. 


1987 


Beuther et al. 


2008). This is a curious chemical differen ce, given the prox i mity of these 


two regions, which has motivated several theoretical studies (IMillar et al. 


1991 


Caselli et al. 


1993 


Rodgers & Charnley 


2001). High spatial resolution maps, which reveal substantial overlap 


between O and N-bearing organics (though the emission peaks are distinct), complicate this 


picture by showing the 

)oundaries between O and b 

-bearing organic rich gas are not well defined 

(Friedel & Snyder 

2008 

Friedel & Widicus Weaver 

2012 


Widicus Weaver & Friedel 

2012 

). This 


'Herschel is an ESA space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from NASA. 
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suggests that differences in the physical conditions (i.e. temperature, density, and gas kinematics) 
on small spatial scales play an important role in the complex morphology of organics in this region 


(see e.g. 


Friedel & Widicus Weaver 


20121) . Additionally, Orion KL harbors at least one outflow, 


often referred to as the plateau, which contains the complex organic CH 3 CN (Wang et ah 


Bell et al. 


2010b 


2014). A detailed picture of the physical conditions from which each molecule emits 


toward Orion KL and other high-mass star forming regions, therefore, provides critical constraints 
for theoretical studies that model the chemistry of complex organics in these environments. 


Unbiased spectral line surveys are useful tools for obtaining the chemica 


inventory of 


Crockett et al 


interstellar molecular gas and constraining the excitation of detected species. 

(2014b, hereafter Paper I) presents a global analysis of the Herschel /HIFI 1.2 THz wide spectral 


survey of Orion KL. These observations were taken as part of the guaranteed time program entitled 
Herschel Observations of EXtra Ordinary Sources (HEXOS). Additional details concerning this 


program can be found in 


Bergin et al. 


( 2010t) . As part of the analysis, Paper I presents models, 


which assume local thermodynamic equilibrium (LTE), of the molecular emission detected in the 
HIFI survey. These models provide estimates for the rotation temperature, T rot , and total column 
density, A [ot , (from which abundances are computed) of each species. In this follow up study, we 
discuss the chemical implications of the LTE models focusing on complex organics. 


This paper is organized as follows. In Sectional we review the Orion KL HIFI scan and LTE 
models. In Section [3] we identify the molecules emitting in the hottest gas using a novel method 
which employs the emission predicted by our spectral models as a direct metric (13.11) and present 
gas-grain chemical models to help interpret our results (13.21) . We discuss our results in Section |4] 
and present our conclusions in Section [5) 
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2. Observations and The Full Band Model 


We obtained spectral scans in all available HIFI (de Gr aauw et al. 20101) bands (la - 7b) 
between 1 March 2010 and 19 March 2011 using the wide band spectrometer (WBS). The full 
spectral survey thus covers a frequency range from 480 to 1907 GHz with two gaps at 1280 - 
1430 GHz and 1540 - 1570 GHz. The WBS provides a uniform spectral resolution of 1.1 MHz 
(0.2 - 0.7 km s" 1 ); the different spatial/velocity components associated with Orion KL are 
thus easily resolved by HIFI. All of the observations were taken in dual beam switch (DBS) 
mode. At HIFI frequencies, the beam size, 6b, of Herschel varies between 11" and 44". In 
bands 1 - 5 (fib = 17"- 44") the telescope was pointed toward coordinates a 72000 = 5 /? 35 m 14.3 s , 

8 72000 = -5°22'33.7", midway between the hot core and compact ridge. In bands 6 and 7, because 
of the smaller beam size (6b = 11" - 15"), separate observations were obtained toward the hot 
core (072000 = 5*35'" 14.5 s , d /2000 = -5°22'30.9") and compact ridge (a mm = 5 /! 35" ! 14.1 s , 


d /2000 = -5°22'36.5"). The data were reduced using the standard HIPE (Ott 


2010 ) pipeline 


(version 5.0, build 1648), which produces “Level 2” flux calibrated double side band (DSB) 
spectra. HIFI provides separate horizontal (H) and vertical (V) polarizations. The DSB spectra for 


both polarizations were c 


sideband (SSB) spectra (Comito & Sch ilke 


econvolved using the doDeconvolution task within HIPE to create single 


20021) after the DSB scans were baseline subtracted 


and spurious spectral features, not flagged by the pipeline, were removed. As a final step, we 
averaged the H and V polarizations together. 

We modeled the molecular emission using the XCLASfj 2 ! program. In total, we detect 
emission and/or absorption from 79 isotopologues originating from 39 molecules, most of which 
have their own molecular fits (those species observed predominantly in absorption or with only 
a few detected transitions over a limited range in excitation energy were not fit). Summing the 
predicted emission from all detected species yields the “full band model”. The fits are constrained 


2 http://www.astro.uni-koeln.de/projects/schilke/XCLASS 
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at mm wavelengths by combining the HIFI spectrum with a ground-based line survey obtained 


with the IRAM 30 m telescope (Terccro et al 


20100 . This survey spans frequency ranges of 80 


- 116 GHz, 130 - 180 GHz, and 200 - 280 GHz, corresponding to atmospheric windows at 
3 mm, 2 mm, and 1.3 mm, respectively. The models are thus constrained from the mm to the 
far-IR and, in almost all cases, include transitions at or close to ground up to excitation energies 
where emission is no longer detected. A reduced^ 2 analysis shows generally excellent agreement 
between the data and models. This is particularly true for complex organics, which typically have 
reduced^ 2 values <1.5 and were constrained by hundreds of lines detected at or above the 3 cr 
level. 


The full Orion KL HIFI spectral survey along with a more detailed description of the data 


reduction procedure and modeling methodology are presented in Paper I and 


Crockett et al. 


(I2014at) . In this study, we use the molecular fits as templates for the data, analyzing the emission 


from the models rather than the data. Examining the HIFI scan in this way has two advantages. 
First, we are able to easily separate emission from different spatial/velocity components. And 
second, we do not have to be concerned with line blends because we examine the model emission 
on a per molecule basis. 


3. Results 

3.1. Characterizing the Molecular Emission 

We employ a novel diagnostic to characterize the molecular line emission within the 
Orion KF HIFI scan. For each molecule and spatial/velocity component, we sum the integrated 
intensity from all T ml sub-components presented in Paper I. We also include emission from 
vibrationally excited modes if they are detected for a given species. We then compute the fraction 
of total integrated intensity originating from states over different ranges of upper state energy, 
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E up . Figs. [TJ [2l and [3] plot these fractions as pie charts for molecules detected toward the hot core, 
compact ridge, and plateau, respectively. We do not discuss emission from the extended ridge 
here because no complex organics are detected toward this component. Molecule names and 
rotation temperatures derived in Paper I are given above each pie chart. If a pie chart includes a 
vibrationally excited mode, it is also indicated in the molecule label. The charts are organized 
so that r rot decreases from the upper left to the lower right. In instances when more than one 
isotopologue is detected for a given molecule, we represent the emission using a rarer species 
to avoid issues with high optical depth. However, when vibrationally excited modes are only 
detected in more abundant isotopologues (e.g. CH 3 CN, S0 2 , etc.), we adopt the more abundant 
species so that emission from vibrationally excited states can be included in the pie charts. For the 
sake of simplicity, we use only main isotopologue names in the discussion below even if the pie 
chart for a particular molecule is made using the spectrum of a rarer species. 

The above metric is a more comprehensive diagnostic for the excitation of a given molecule 
compared to a single T mt because many species had to be fit with multiple temperature sub¬ 
components indicating the presence of temperature gradients. Such molecules therefore probe a 
range of temperature environments within a given spatial/velocity component. Moreover, in order 
to accurately fit the vibrationally excited emission of some species, we had to adjust T rot , and, in 
some instances, increase jV t ot relative to the ground vibrational state, suggesting these molecules 
are being pumped from mid/far-IR photons. 

Examining Figs. |T] — [3] we see that those molecules emitting a significant fraction of their 
total line radiation from high energy states (.E U p > 800 K) are readily apparent. Such molecules 
represent the most highly excited molecules within Orion KL and likely originate from the 
hottest gas. We emphasize that this metric is independent of rotation temperature and relies 
solely on the observed emission as predicted by our models. Ethyl cyanide (C2H5CN) is the only 
cyanide detected toward the hot core with a rotation temperature less than 200 K (7j. ot = 136 K). 


- 8 - 


Its LTE model, however, requires a hotter, more compact sub-component (T mt = 300 K) to 
accurately fit the emission (Paper I). The presence of this sub-component is evidenced in the 
pie chart of C 2 H 5 CN because approximately 5% of its integrated intensity originates from states 
with E U p > 800 K. This is in contrast to molecules with similar rotation temperatures (e.g. 
H 2 CO), which do not emit significantly from states with E up > 800 K. Methyl cyanide (CH 3 CN), 
furthermore, does not have an especially high r rot toward the plateau suggesting it does not probe 
particularly hot gas in the outflow. Including emission from the v 8 = 1 vibrationally excited mode 
(Eup > 520 K), however, indicates that it may be probing hotter material than suggested by the 
derived T ml . 

Of the 39 molecules detected in the HIFI spectrum, 8 are complex. This subset includes 
species which contain only carbon, hydrogen, and oxygen as well as nitrogen bearing organics 
which may also contain oxygen (i.e. NH 2 CHO). From Fig. [[] we see that among the 5 complex 
molecules detected in the hot core, the nitrogen bearing species CH 3 CN, C 2 H 3 CN, and C 2 H 5 CN, 
all of which are cyanides, display more highly excited emission, as represented by the pie charts, 
than their oxygen bearing counterparts, CH 3 OH and CH 3 OCH 3 . Moreover, the simpler cyanides 
HCN, HNC, and HC 3 N also display some of the most highly excited emission toward the hot 
core. Fig. [2] shows a similar trend in the compact ridge. The N-bearing complex organics CH 3 CN 
and NFFCHO display more highly excited emission than organics which contain only C, H, and 
O: CH 3 OH, C 2 H 5 OH, CH 3 OCH 3 , and CH 3 OCHO. Finally, Fig. [3] shows that CH 3 CN, the only 
complex organic detected toward the plateau, displays the most highly excited spectrum. 

Figs. |Tj-[3] also show that sulfur bearing molecules which also contain oxygen, i.e. SO, S0 2 , 
and OCS, make up another subset of molecules that trace the hottest gas toward the hot core, 
compact ridge, and plateau. The two exceptions to this are S0 2 detected toward the compact 
ridge (T ml = 100 K) and OCS detected toward the plateau (T rot =110 K). We note, the rotation 
temperature of S0 2 toward the compact ridge is not well constrained because the line profiles are 
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dominated by the plateau and hot core and thus could be underestimated. The relatively low T ml 
value for OCS, combined with its pie chart, however, suggest it is probing cooler material than SO 
and S0 2 toward the plateau. 


3.2. Chemical Modeling 


In order to determine if theoretical studies agree with our derived abundances and T Iot values, 


(Garrod & Herbst 


2006; 

Garrod et al. 

2007, 

OO 

o 

o 

<N 


exist which include 1-D hydrodynamical treatments (e.g. 


Aikawa et al. 


2 , 


20121) . These models, 


however, are not easily compared with our single dish results because the predicted molecular 
abundances depend upon a time dependent physical profile (i.e. temperature, density, and velocity 
field as a function of radius). Because we do not attempt to model the physical structure or 
dynamical history of Orion KL here, we use the simpler “warm-up” approach described below. 
We emphasize that a global chemical model that reproduces the abundances of all observed 
species in the Orion KL HIFI survey is beyond the scope of this work because predicted gas 
phase abundances are sensitive to the details of the thermal history for a given region, binding and 
desorption energies of grain surface molecules, and reaction branching ratios, all of which can be 
highly uncertain. Rather, we aim here to determine if the OSU gas-grain chemical code, which 
assumes a relatively simple thermal history (see below), predicts abundances for complex organics 
that are comparable to what is observed in the HIFI scan and if this model can reproduce the 
observed trend of N-bearing complex molecules probing hotter gas relative to O-bearing organics. 


The OSU gas-grain chemical model assumes a two phase evolution. The first phase is 
an isothermal collapse in which the gas density increases from 3 x 10 3 to 1 x 10 7 cm~ 3 at 
a temperature of 10 K over ~10 6 years. The second is a “warm-up” phase in which the gas 
temperature increases by following a quadratic power law over a warm-up timescale, t warm , to 
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a final temperature, T F . The gas and dust are assumed to be well coupled to one another, thus 
the temperatures of both components are assumed to be equal. Within this model, complex 
organics are formed both in the gas phase and on grain surfaces. Many of the complex organics 
are produced by reactions between radicals, which are formed primarily by photodissociating 
molecules abundant in ice mantles, during the warm-up phase. Because the warm-up stage is 
characterized by high dust extinction, most of the photodissociation is the result of photons 


induced by cosmic rays (Ruffle & Herbst 


2001). The smooth rise in temperature over time is 


therefore key because it allows heavy radicals to become mobile on grain surfaces without being 
sublimated, and facilitates gas phase pathways with activation energies. For the grid, we adopt 
two values for T F of 200 K and 300 K, commensurate with the highest T mt values measured in the 


and intermediate mass star formation as laid out by 

Viti et al. 

( 

2004) 

. All other model parameters, 

including the chemical network, are the same as 

Garrod et 

al. 

(200* 

1). We also adopt the same 

initial fractional abundances from the “standard” 

model of 

Garrod et al. 

(2008 

). 


Figures 0] and [5] plot the predicted gas phase abundances of the eight complex molecules 
detected in the HIFI survey as a function of time during the warm-up phase for t warm = 5 x 10 4 years 
and 2 x 10 5 years, respectively. Measured abundances for these species toward the hot core and 
compact ridge (Paper I) are also overlaid as horizontal lines. Table |T| reports the temperature 
at which the abundance of each molecule is maximized, T max , along with the peak abundance. 
Examining Table |U we see that changing 7) or / warm does not strongly affect T max , which varies 
by < 11% except when a molecule’s peak abundance occurs at or just before T F , in which 
case T max « T v . The gas-phase abundances of those species with T max < T F typically peak at 
temperatures somewhat higher than their sublimation temperature when they thermally desorb 
from grains, which is why T max does not vary much for these molecules. The efficient sublimation 
of a molecule is shown in Figures 0] and [5] as a steep rise in gas phase abundance followed by a 
slow decrease corresponding to destruction pathways which are active at high temperatures. We 
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ref er here to rises in abundance that occur last in the evolution of the models, but note additional 
earlier abundance peaks occur for many of the organics which are the result of interdependent gas 
phase and grain surface pathways. Those species with T max « T v ., on the other hand, (i.e. C 2 H 3 CN, 
NH 2 CHO, and CH3OCH3) have formation/destruction pathways which yield net increases in their 
abundances after sublimation from grain surfaces. The peak abundance varies more significantly 
than T max across the grid (by as much as a factor of ~ 7), especially when t warm is changed. As 
noted previously, the OSU code typically predicts higher abundances for complex organics when 
/warm is longer because there is more t ime at intermediate temperatures (~ 20 - 40 K) when heavy 


2006; 


Garrod et al 


2008). The parameter 


radicals are mobile on grain surfaces (Garrod & Herbst! 
t w thus plays the more important role relative to T v , in the production of complex molecules. 


We include in our analysis the N-bearing organics C 2 H 3 CN and C 2 H 5 CN which have not 
been discussed in earlier studies that employ the OSU chemical network. Within the model, 
C 2 H 3 CN forms predominantly in the gas phase via the reaction, 

CN + C 2 H 4 —> C 2 H 3 CN + H, (1) 

though dissociative recombination of C 3 H 4 N + contributes somewhat when temperatures are 
< 40 K. The main destruction routes are through the ion-neutral reaction, 

HCO + + C 2 H 3 CN -> C 3 H 4 N + + CO, (2) 

and accretion onto dust grains. At temperatures > 40 K accretion onto dust grains is the major 
destruction mechanism for gas-phase C 2 H 3 CN, after which C 2 H 3 CN is hydrogenated on grain 
surfaces (7\ in < 120 K) or quickly thermally desorbed back into the gas phase (7) (m > 120 K). 
The formation of C 2 H 5 CN occurs on grain surfaces via multiple hydrogenations of HC 3 N 
and C 2 H 3 CN. C 2 H 5 CN remains on grain surfaces until it is efficiently thermally sublimated at 
temperatures > 100 K, at which point the gas phase abundance rises rapidly. 


Comparing the model predictions with observations, we see that the chemical models 
reproduce observed abundance levels for a majority of the complex organics at some point in 
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time within the grid. The species C 2 H 3 CN, C 2 H 5 CN, and CH3OCHO, however, are consistently 
under predicted, though the agreement is better for the longer / wami , which produce systematically 
higher abundances for these species than the shorter t waim • We note that the abundances for these 
molecules are under predicted by less than 3cr for f warm = 2 x 10 5 years given the uncertainties 
in the observations (Paper I). Moreover, abundances for CH3CN, toward the hot core and 
compact ridge, as well as CH 3 OCH 3 , toward the compact ridge, are only under predicted for 
hvarm = 5 x 10 4 years. Hence, we conclude models corresponding to t warm = 2 x 10 5 years, in 
general, produce abundances that agree more closely with the observed values. For these models, 
observed abundances are reached at times > 10 5 years. We note that our model grid predicts peak 
abundances for NH 2 CHO that are more than three orders of magnitude higher than that observed 
toward the compact ridge and even approximately 2 orders of magnitude higher than the NH?CHO 


abundance measured toward SagB2(N) (Neill et ah 


20140 . This overabundance is a result of the 


fact that NH 2 CHO forms efficiently in the gas phase by the combination of NH 2 and H 2 CO, 
and is continually accreted onto and desorbed from dust grains without an efficient destruction 
mechanism late in the model’s evolution. Currently the major gas-phase destruction mechanism 
(behind accretion onto dust grains) is via reactions with H 3 0 + and H ! that form NH 2 CH 2 0 + , 
which is subsequently destroyed by dissociative recombination producing NH 2 , HCO, and H 2 CO. 
Although the gas phase route helps to explain the presence of NH 2 CHO in hot gas, we emphasize 
that its efficiency is quite uncertain and likely too large to explain abundances observed toward 
hot cores. 


4. Discussion 
4.1. Chemical Implications 

Our results indicate that complex N-bearing organics, i.e. cyanides and NH 2 CHO, probe 
hotter environments than complex species which contain no nitrogen. Previous single dish and 
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interferometric studies have show n that gas which traces the hot core is hotter than gas which 


traces the compact ridge (see e.g. 


Blake et al. 


1987 


Beuther et al. 


20051). Because the hot core 


is richer in N-bearing organics, these molecules tend to trace hotter gas than their O-bearing 


counterparts. Recently 


Friedel & Widicus Weave r (2012) overlaid high spatial resolution maps 


of key N/O-bearing complex molecules with a temperature map constructed from a combination 


of their own CH 3 OH data and Nth observations from 


Goddi et al. 


(2011). They showed that 


the N-bearing species C 2 H 5 CN, C 2 H 3 CN, and CH 3 CN are cospatial with one another and trace 
higher temperature regions associated with the hot core and nearby sources SMA1 and IRc7. 


The O-bearing mo 


Favre et al. 


201 lbi: 


ecules CH 3 OCHO and CH 3 OCH 3 , also cospatial with one another (see also 


Brouillet et al. 


20131) . have emission peaks near cooler regions associated with 


the compact ridge and sources IRc5 and IRc6. We emphasize that our results apply to molecules 
detected within a spatial/velocity component, i.e. N-bearing organics within the hot core are hotter 
relative to O-bearing molecules within the hot core with the same being true for the compact ridge. 
Our excitation temperatures, furthermore, are derived individually from the molecular emission 
of each species, and are the most robust to date because the emission is constrained by hundreds 
to thousands of transitions over mm, sub-mm, and far-IR wavelengths. The inclusion of up to 
thousands of lines per molecule in our analysis thus places these results on a strong statistical 
footing. 

If both O and N-bearing complex organics form predominantly on grain surfaces, cyanides 
along with NH 2 CHO may be more difficult to remove from grain surfaces than O-bearing 
species. Within a hot core, we might then expect oxygen bearing organics to be released 
during an earlier, presumably cooler, epoch and/or further from the central protostar. Measured 
excitation temperatures should thus be higher toward the same location, and spatial temperature 
variations larger and clumpier for N-bearing organics relative to O-bearing species because 
the former traces hotter material along the line of sight. Because interferometric observations 
show Orion KL has structure on spatial scales < 1" (e.g. 


Wang et al. 


2010 ; 


Goddi et al 


2011 
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Favre et al.11201 lac LFriedel & Widicus Weaverll20111) . our results illustrate the need for excitation 
temperature maps derived from both O and N-bearing organics at sub-arcsecond resolution 
to see if differences exist. Such observations will surely be attainable with ALMA. Existing 


interfero metric maps of the kinetic tempe rature derived from CH 3 CN (IWang et al. 


CH 3 OH (Friedel & Widicus Weavei 


20101 ) and 


2012 ) tentatively support this viewpoint at least toward the 


compact ridge where r rot values are measured to be ~ 170 - 280 K for CH 3 CN and < 200 K for 
CH 3 OH. We, however, note that the spatial resolution of the CH 3 CN observations (1.6" x 0.9") is 
about a factor of two higher relative to the CH 3 OH map (2.2" x 2.0"). 

Another possibility is that hot gas phase chemistry may be producing the highly excited 


cyanides in the hot core. 


Wang et al. 


(120101 ) . for example, argue that Cf 


forms in the gas based on the chemical models of 


CN in the hot core likely 


Rodgers & Cham lev (2001). These models 


show that at high temperatures (7 kin = 300 K), hot cores produce CH 3 CN, HC 3 N and HCN in the 
gas phase at observed abundance levels on timescales of < 10 5 - 10 6 years. Models presented by 


Caselli et al. 


(119931) show similar results. In this scenario, cyanides naturally trace hotter materia l 


because they form efficiently in the gas phase at high temperatures. 


Hemandez-Hemandez et al. 


(2014), furthermore, presents a correlation in the abundance of CH 3 CN with temperature based 


on LTE modeling of the 12^-11^ ladder toward 17 hot cores, which suggests hotter material 
facilitates methyl cyanide formation, presumably in the gas phase, because measured temperatures 
in the hottest component of each model (2 temperature components were used for each core) range 


from approximately 100 - 500 K. 


Gerin et al. 


(1992) reports a CH 2 DCN / CH 3 CN ratio > 0.005 


towar d the Orion K L hot core, comm ensurate with D/H ratios measured for FRO (INeill et al. 


2013b) and CH 3 OH (Neill et al. 


2013aj) . suggesting at least some methyl cyanide forms via grain 


surface chemistry at low temperatures. If gas phase formation routes are active for N-bearing 
organics in the hottest gas, we would expect gradients in the D/H ratios of complex N-bearing 
species like methyl cyanide. 
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Previous observational studies have investigated possible differentiation between O and 
N-bearing organics toward other hot cores. Using single dish spectroscopic observations in the 


mm and sub-mm, 


Bisschop et al.l (2 0071) computed the abundance and rotation temperatures 


for several O and N-bearing organics toward a sample of 7 hot cores. They found that the 
abundances between the complex “hot”, which they define as T mt > 100 K, O-bearing species 
CH 3 OCH 3 , CH 3 OH, C 2 H 5 OH, CH 3 OCHO, and H 2 CO are correlated with one another. However, 
abundances between the O-bearing and the N-bearing molecules HNCO and NH 2 CHO, both of 
which trace hotter gas than the O-rich organics toward Orion KL (Fig. [Hand |2] respectively), 


did not correlate. 


Bisschop et al. 


(2 0070 argue that these results indicate that complex O and 


N-bearing organics originate from grain surface chemistries that are not strongly coupled to 


one another. Conversely, 


Fontani et al.1 ([ 20071) carried out a similar study focusing on C 2 H 5 CN, 


C 2 H 3 CN, and CH3OCH3 using a sample of 6 hot cores. They found that abundances for all 
three species correlated with one another indicating that a simple oxygen/nitrogen chemical 
dichotomy is likely an over simplification. Indeed several interferometric studies have noted 
that (CH 3 ) 2 CO has a unique morphology toward Orion KL, peaking primarily toward regions 


where emission from both O and N-bearing complex organics are detect ed (IFriedel & Snvde 


2008 


Friedel & Widicus Weaver 


20121 : 


Widicus Weaver & Friedel 


2012). C 2 H 5 OH also has a 


unique morphology toward Orion KL peaking at or near the (CH 3 ) 2 CO peak as opposed to tracing 


CH3OCHO and CH3OCH3 (Guelin et al 


2008). Many of these studies have also noted that there 


is significant overlap in the spatial distributions of large organics that contain O and N, indicating 
that chemically differentiation toward Orion KL, in the strictest sense, is too elementary a picture. 

Our results also show that S + O-bearing molecules probe environments as hot as N-bearing 
complex organics, suggesting the formation pathways for these species favor high densities 
and temperatures. Most theoretical studies which trace sulf ur chemistry in hot cores assume 


H?S forms on 


Hatchell et al. 


grain surfaces via hydrogenation of S atoms (Miliar et ah 


1997 


Charnlev 


19971 


1998). Within these models, all other S-bearing molecules are produced in the gas 
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phase after H?S is evaporated from grain surfaces. Once in the gas, H?S is efficiently destroyed by 
H atoms or other radicals or ions such as H 3 0 + with most of the sulfur eventually channeled into 
SO and S0 2 through reactions with OH and atomic O. Paper I reports H 2 S hot core abundances 
that are higher than that of SO and SO? by factors of 20 and 8, respectively, suggesting H?S is 
being actively replenished via evaporation from icy grain mantles. The high rotation temperatures 
of SO, SO?, and OCS, thus, may mean that sulfur is being most efficiently channeled from H?S 
into these species within the hottest gas. The fact that OCS traces cooler gas than SO and SO? 
toward the plateau, however, may indicate a unique formation route within the outflow for this 
species, possibly related to shock chemistry. 


4.2. Comparison with Chemical Models 


Do the models presented in Section 13.21 which assume a rather simple warm-up phase, 
predict that N-bearing organics should systematically probe hotter gas than O-bearing organics? 


We take the same approach as 


T max (Tabled]). As noted in 


Garrod et al 


Garrod et al. 


(2008) and compare the observed T rol to the predicted 


( 20081) . the OSU chemical code computes abundances for 


a single point in space over time. However, molecules with modeled abundances that peak at high 
temperatures represent better tracers of hot gas. Orion KL is a complex region with temperature 
and density gradients. Different species probe different environments, with the hottest and densest 
gas corresponding to material in a more evolved stage in its warm-up progression relative to 
cooler, more extended material. The comparison of T ml with T max is thus an efficient way to 
determine if the chemical models predict N-bearing organics are tracers of the hottest gas. 


Figure [6] plots T rot for the complex organics detected in the HIFI survey with T max values 
overlaid as a function of molecule toward the hot core and compact ridge. For C?H 5 CN, we plot 
the r rot value of the compact 300 K component required to fit the emission (see Sec. 13.11) as a 
green square to indicate that this species probes gas as hot as the other cyanides toward the hot 
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core. We note that abundances plotted in Figures |4] and [5] for most species stay relatively flat late 
in their evolution, which yield curves that do not have sharp peaks. Consequently, these molecules 
are predicted to be abundant over a range of temperatures, but are most abundant either right 
after sublimation (i.e. CH 3 CN, C2H5CN, CH3OH, C9H5OH, and CH3OCHO) or continue to rise 
in abundance as the gas gets hotter (i.e. C2H3CN, NH 2 CHO, and CH3OCH3). Our values for 
T max are thus not meant to be robust predictions for T m . However, because the high molecular 
abundances observed toward Orion KL are primarily due to active grain mantle sublimation, we 
use T max as a metric to, first, identify which molecules are predicted to be most abundant in the 
hottest gas, second, determine if there is general agreement between the observations and our 
models, and, third, identify the molecules that are obviously discrepant. 

Figure [6| shows that most species have T max / T ml combinations that agree to within 
2cr. However, CH 3 CN and CH3OCH3 have T max predictions that are markedly different from 
their corresponding T vot values. For CH 3 CN, we observe T rol > 230 K, but the models yield 
Tmax < 111 K. As noted above, however, the peak CH3CN abundance, which occurs shortly after 
it thermally desorbs from grain surfaces, is not very sharp, and there is almost as much gaseous 
CH 3 CN predicted at temperatures of ~ 200 K - 300 K. For CH3OCH3, we observe T mt ~ 100 - 
110 K, but the models predict T mt > 200 K. This is a result of the fact that CH3OCH3 continues to 
be produced in the gas phase, after thermal desorption at approximately 70 K, via the reaction, 

CH3OH+ + CH3OH CH3OCH+ + H 2 0 (3) 

and subsequent dissociative recombination. The formation of CH3OCH3 is thus sensitive to 
the gas phase abundance of methanol which is over predicted in our models at its peak. The 
contribution of Reaction [3] to the CH3OCH3 abundance may, therefore, be over estimated in our 
models. We note that the predicted T max for C 2 H 5 CN agrees well with our measured 7 j. ot , but is 
much lower than 300 K, the temperature of the hot, compact component required to fit the HIFI 
emission (see Sec. 13.11) . which is plotted as a green square in Figure [6l However, just as with 
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CH 3 CN, the peak abundance for C 2 H 5 CN is not sharp and similar abundances are predicted for 
temperatures > 125 K. 


The age of the hot core is estimated to be approximately 10 3 -10 4 years based on measu red 


(Plambeck & Wright 

1987; 

Walmslev et al. 

'O 

00 

Brown et al. 


1988). The models presented in this study that have f warm = 2 x 10 5 years, which produce 


better agreement overall with measured abundances of complex organics relative to models with 
t Wia rm = 5 x 10 4 years, require timescales > 10 5 years to reproduce observed abundance levels (see 
Section lT2l) . The time necessary to form complex organics in sufficient quantities according to 


our models is thus at least an order of magnitude larger than t 


by D/H ratios of lighter hydrides. Moreover, 


le age of the hot core as estimated 


Esplugues et al 


(2014), who model the sulfur 
chemistry toward the Orion KL hot core in detail, require an intermediate age of ~5 x 10 4 years to 
reproduce the observed abundances of S-bearing molecules in their best fitting models. We note 
that these timescales are not well constrained given the complex density/temperature structure 
and the uncertain dynamical history of Orion KL. As a result, we do not consider our timescale 
of > 10 5 years as a robust age estimate for the hot core or compact ridge. The fact that it is 
significantly longer than other hot core age estimates, however, suggests certain events, not 
modeled here, may have facilitated the production of complex molecules toward Orion KL early 
in its evolution. Lor example, if the gas in the hot core and compact ridge spent a significant 
amount of time at temperatures ~ 20 - 40 K prior to warmup, the formation of complex organics 
may have been hastened due to the mobility of heavy radicals on grain surfaces at earlier times. 
Such a possibility is not unreasonable given that most molecules measured toward the extended 
ridge have measured rotation temperatures > 20 K (Paper I). The clumpy structure of Orion KL 
may also have accelerated the formation of complex organics because holes in dusty material 
could have provided pathways for UV photons to break up the molecular constituents of ice 
mantles, facilitating the formation of heavy radicals on dust grains. These caveats may also help 
to reconcile the under predictions of C 2 H 3 CN, C 2 H 5 CN, and CH 3 OCHO by our model grid. We 
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also note that if the 
than 720 years ago 


lot core is indeed heated externally by an explosive event which occurred 


Bally 


2008 


Zapata et al. 


2009. 


2011 


Bally et al. 


2011 


Nissen et al. 


2012 ), 


ess 


then the dynamical history of Orion KL is far more complex than simulated here, which may also 
help to explain why some discrepancies exist between measured and predicted abundance levels. 


5. Conclusions 

Using spectral models fit to the molecular emission within the Orion KL HerschelfHlFl 
spectral survey, we have identified the complex organics that are emitting in the hottest gas 
toward the different spatial/velocity components associated with this region. We find that the 
nitrogen bearing complex organics CH 3 CN, C 2 H 3 CN, C 2 H 5 CN, and NH 2 CHO systematically 
trace hotter gas than O-bearing complex species which contain no nitrogen: CH 3 OH, C 2 H 5 OH, 
CH 3 OCH 3 , and CH 3 OCHO. We emphasize that these results apply to molecules detected within a 
spatial/velocity component and are obtained using the emission predicted by our models, which 
reproduce the data well, as a direct diagnostic of the excitation. Hence, they do not depend on any 
model assumptions (e.g. LTE). If complex organics form on grain surfaces, these observations 
may indicate that complex N-bearing molecules are more difficult to remove from grain surfaces 
than organics which contain only H, C, and O. Another possibility is that the complex N-bearing 
species are forming via a high temperature gas phase chemistry and thus naturally trace hotter gas. 

We also compared our derived molecular abundances and rotation temperatures to chemical 
models derived using the OSU gas-grain chemical code, which includes both an isothermal 
collapse and warm-up phase. We find better agreement in the observed abundance levels of 
complex organics when / wami is set to 2 x 10 5 years as opposed to 5 x 10 4 years. Abundances 
for C 2 H 3 CN, C 2 H 5 CN, and CH 3 OCHO are under predicted for all models in our grid, though 
the model and observed abundances disagree by less than 3cr for / warm = 2 x 10 5 years. We 
also find that model 7 niax values agree reasonably well (within 2cr) with observed rotation 
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temperatures for all complex organics except CH 3 CN and CH3OCH3. Finally, SO, S 0 2 , and 
OCS, i.e. S + O-bearing molecules, are among the molecules which trace the hottest gas toward 
Orion KL. If these species are formed via gas phase reactions starting with the sublimation of H 2 S 
from grains, our results may indicate that the formation routes of these species are most active at 
high temperatures. 
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Table 1. Model Predictions for Complex Molecules 


Tf 

^warm 

200 K 

5 x 10 4 yr 

200 K 

2 x 10 5 yr 

300 K 

5 x 10 4 yr 

300 K 

2 x 10 5 yr 

Molecule 

«x/«K 2 

T m ax (K) 

n x /n. h. 

T m ax (K) 

n x /na 2 

T m ax (K) 

Ux/«H 2 

T m ax (K) 

CH3CN 

7.9 x 10~ 9 

104 

3.5 x 10' 8 

101 

6.8 x 10~ 9 

111 

2.8 x 10~ 8 

100 

C 2 H 3 CN 

1.9 x 10“ 9 

200 

3.4 x 10- 9 

200 

2.3 x 10~ 9 

300 

3.2 x 10~ 9 

300 

C2H5CN 

8.4 x 10~ 9 

127 

4.6 x 10~ 8 

123 

7.0 x 10~ 9 

126 

2.9 x 10 -8 

122 

NH 2 CHO 

2.3 x 10" 6 

200 

4.1 x 10~ 6 

200 

3.1 x 10~ 6 

300 

5.0 x 10" 6 

295 

CH 3 OH 

4.8 x 10~ 5 

119 

4.0 x 10~ 5 

115 

4.9 x 10- 5 

118 

4.1 x 10~ 5 

122 

C2H5OH 

5.0 x 10~ 8 

127 

1.8 x 10- 7 

123 

3.8 x 10 -8 

126 

1.5 x 10- 7 

122 

CH 3 OCH 3 

9.3 x 10~ 8 

200 

2.5 x 10- 7 

200 

1.1 x 10- 7 

300 

2.9 x 10- 7 

300 

ch 3 ocho 

6.6 x 10~ 9 

92 

1.1 x 10~ 8 

83 

6.8 x 10- 9 

91 

9.5 x 10~ 9 

88 


Note. — Abundances represent peak values. 
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Fig. 1.— Pie charts that plot the fraction of total integrated intensity originating from states in E up 
ranges 0 - 200 K (dark blue), 200 - 800 K (cyan), and 800 - 3000 K (red) for molecules detected 
toward the hot core. The molecule ID and T mt derived in Paper I are given above each chart. If 
a chart includes emission from a vibrationally excited state, it is indicated in parentheses with the 
molecule ID. 
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Fig. 2.— Same as Fig.|T|for the compact ridge. 
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Fig. 3.— Same as Fig.|T|for the plateau. 
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Fig. 4.— Predicted gas phase abundances with respect to FF for the eight complex organics de¬ 
tected in the Orion KL HIFI survey plotted as a function of time for t warm = 5 x 10 4 years. The 
abundance of each species is plotted as a different color (solid lines). Abundances for N and O- 
bearing organics are plotted in the top and bottom rows, and models with 7) = 200 and 300 K are 
plotted in the left and right columns, both respectively. Observed abundances toward the hot core 
and compact ridge are plotted as horizontal dashed and dotted lines, respectively, with the same 
color scheme as the model values. The bottom x-axis is time such that zero corresponds to the 
beginning of the warm-up phase and the top x-axis displays the corresponding temperature. 
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Fig. 5. 


Same as Figure|4]for / warm = 2 x 10 5 years. 
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Fig. 6.— Rotation temperatures (black dots) measured toward the hot core (left panel) and compact 
ridge (right panel) plotted as a function of molecule with T max values from our model grid (Tabled)) 
overlaid. The range of T max predicted over all grid models for each molecule is represented by a 
red bar if T max < T v , and a red triangle at 200 K if T max w 7) indicating these molecules probe 
temperatures > 200 K in our model grid (see Section 13.21 ). The green square represents the hot, 
compact component required to fit the C 2 H 5 CN emission (see Section I3T1) in the HIFI scan. The 
error bars on the observed r rot values represent ler uncertainties. 
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